Long-range steady state density profiles induced by localized drive 
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We show that the presence of a localized drive in an otherwise diffusive system results in steady- 
state density and current profiles that decay algebraically to their global average value, away from 
the drive in two or higher dimensions. An analogy to an electrostatic problem is established, 
whereby the density profile induced by a driving bond maps onto the electrostatic potential due 
to an electric dipole located along the bond. The dipole strength is proportional to the drive, 
and is determined self-consistently by solving the electrostatic problem. The profile resulting from 
a localized configuration of more than one driving bond can be straightforwardly determined by 
the superposition principle of electrostatics. This picture is shown to hold even in the presence of 
exclusion interaction between particles. 

PACS numbers: 05.40.-a, 05.70.Ln, 05.40.Fb 



in 



i 

C 

O 

o 



> 
00 

m 
oo 

O 



X 



The effect of a local perturbation on the steady state 
density profile of systems of interacting particles has been 
studied in a wide variety of contexts. When a system 
under thermal equilibrium conditions is perturbed by a 
localized external potential, the equilibrium density is 
changed only locally under generic conditions. This is a 
result of the fact that as long as the system is not at a 
critical point, it is characterized by a finite correlation 
length. This is not necessarily the case in driven systems 
where detailed balance is not satisfied Algebraically 
decaying correlations in the steady state profiles have 
been found in a number of boundary driven models such 
as the Asymmetric Simple Exclusion Process (ASEP) [2j, 
interface growth models^ and transport models in one 
and higher dimensions [4j-[6j] . 

A natural question is what happens when the drive 
is localized not on the boundaries but in the bulk. In 
a recently studied model of a one-dimensional Symmet- 
ric Simple Exclusion Process (SSEP) it was shown that 
the presence of a single driving bond (battery) in the 
bulk does not generate any algebraic density profile in 
the steady state Q ■ The density profile was found to be 
flat away from the battery, albeit with a discontinuity at 
the location of the battery. 

In the present Letter we consider the effect of driving 
bonds localized in a finite region in an infinitely large 
many body system. We demonstrate that under rather 
generic conditions a localized drive in an otherwise equi- 
librium system in dimensions higher than one, results in 
a steady state density profile with an algebraically decay- 
ing tail. This is done by first studying the case of non- 
interacting particles diffusing on a d-dimensional lattice 
with a directional drive along a single bond (battery). 
We then generalize the results to arbitrary localized con- 
figurations of driving bonds, and to the case of particles 
with exclusion interaction. 

In the case of non-interacting particles with a single 
driving bond, we show that the density profile can be 



mapped onto the electrostatic potential generated by 
an electric dipole located at the driving bond, whose 
strength can be calculated self-consistently. Thus, for 
example, in d = 2 dimensions, the density profile de- 
cays as 1/r at distance r away from the driving bond, in 
all directions except the one perpendicular to the drive, 
where it decays as 1/r 2 . More interestingly, other lo- 
calized configurations of driving bonds result in different 
power-law profiles. In this case the density profile can 
be determined by a linear superposition of the profiles 
generated by each driving bond. For example, when the 
electric dipoles corresponding to two driving bonds form 
a quadrupole, the density profile generically decays as 
1/r 2 while in some specific directions it decays as 1/r 4 , 
at large distances. The correspondence to the electro- 
static problem still holds when local exclusion is switched 
on. The only difference is in the dipole strength which, 
unlike the noninteracting case, can not be determined 
self-consistently. In the interacting case, our results thus 
generalize the one dimensional situation studied in 0] 
and show that for d > 2, the density profile decays alge- 
braically away from the battery. 

We start with the simple case of non-interacting parti- 
cles diffusing in a medium with a single driving bond. 
As an illustration we consider explicitly the d = 2 
case. Generalization to arbitrary dimensions is straight- 
forward. We consider a two-dimensional square lattice 
of L x L sites with periodic boundary conditions. There 
are N = pV noninteracting particles where V — L 2 is 
the number of sites and p is the global average den- 
sity of particles. Each particle performs an indepen- 
dent random walk in continuous time. A particle at site 
f = {m, n) ^ (0, 0) can hop to any of the neighboring 
sites with rate 1. We introduce a localized drive by set- 
ting the hopping rate across the bond (0,0) — > (1,0) to 
be (1 — e) with e < 1 (see FigQJi). In the absence of the 
localized drive (e = 0), detailed balance holds and the 
system reaches a steady state with a flat density profile 
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FIG. 1. (a) Bond-configuration with a single driving bond 
(0,0) — > (1,0) and hopping rate e = (1 — e). Everywhere 
else the rate is 1. The and symbols denote the sign of 
the electric charge. The symbol represents the magnetic 
flux coming out of the plane, and represents the flux going 
in. (b) A configuration of closed loop of driving bonds gen- 
erating no electric charges. The magnetic flux through the 
central plaquette is four times larger than that of the other 
four neighboring plaquettes. 



with density p at each site. When the localized drive 
e is switched on, it manifestly violates the detailed bal- 
ance and induces a global modification of the steady state 
density profile. 

Let 0(m, n, t) denote the average density of particles 
at site (m, n) at time t. Its time evolution can be eas- 
ily written down by counting the incoming and outgoing 
moves from each site. For sites (m, n) ^ {(0,0), (1,0)}, 
it is easy to see 0(m, n, t) satisfies the standard diffusion 
equation, dt<p(in,n,t) = V 2 0(m,n,i), where V 2 is the 
discrete Laplacian: V 2 0(m, n) — <fi(m + l,n) + 0(m — 
1, n) + <f>(m, n + 1) + 0(m, n — 1) — 40(m, n). The evo- 
lution equations are slightly different for the two spe- 
cial sites (0,0) and (1,0) connecting the driving bond: 
d t <t>(0,0,t) = V 2 0(O,O,t) + £0(0,0, t) and 9*0(1, 0,i) = 
V 2 0(1, 0, t) — £0(0, 0, t). It is useful to write these evolu- 
tion equations in a combined form by introducing Kro- 
necker delta symbols. For brevity, let us also denote 
r = (m, n), = (0,0) and ei = (1,0). In the long time 
limit the system approaches a time independent steady 
state 0(f) satisfying 



(f) = -£0(0) 



S r,0 



(i) 



for all f. Here is the Kronecker delta function. It 
is instructive to first note the formal resemblance of ([1]) 
with the Poisson equation (lattice version) of the 2-d elec- 
trostatic problem. One identifies 0(f) as the electrostatic 
potential and the right hand side (rhs) of ([1]) is identi- 
fied with two point charges of equal strength but of op- 
posite signs sitting respectively at the two ends of the 
driving bond. Thus we have effectively a dipole sitting 
on the weak bond. However, unlike in standard electro- 
statics the charge strength £0(0) has to be determined 



self-consistently. 

We consider the thermodynamic limit L — s- oo, ./V — > oo 
with density per site p = N/V fixed. The exact solution 
of ([TJ at any f can be expressed in terms of the lattice 
Green's function G(f, fo) which is the Coulomb potential 
due to a single point charge of unit strength at fo, that 
satisfies V 2 G = —6f ? . Using the superposition principle 
one obtains the solution 



f) =p + £0(O) G(f,0)-G(f,e!) 



(2) 



The constant 0(0) can be determined self-consistently 
by substituting f = in (|2|). This gives, 0(0) = 
p/[\ - £(G(0, 0) - G(0, Si)]. Then by evaluating the lat- 
tice Green's function for an infinite square [8], one finds 
0(O)=p/[l-£/4]. 

To determine the large distance behavior of the so- 
lution, one can use the continuum approximation un- 
der which the Green's function behaves as G(f, fo) ~ 
— 2^1n|f — fo|, for large |f — fo|. Substituting in ©, 
one finds that the density in ([2]) decays for large r alge- 
braically as 



0(f) 



£0(0) ei 
2tt j 



(3) 



This density profile also leads to a nontrivial current pro- 
file. The average particle current density, away from the 
drive, is j(f) = — V0(f) and it decays for large r as 



£0(0) 1 
2tt r 2 



ei 



2(e\ • f)f 



■0(gs)- (4) 



In the electrostatic analogue, j(r) is precisely the electric 
field generated by the dipole. 

Due to the superposition principle, the above analysis 
can be readily generalized to the case of arbitrary local- 
ized configuration of the driving bonds. For example, 
consider a case of two driving bonds with rates (1 — e) 
each, one from (0,0) to (1,0) and the other from (0,0) 
to (—1,0), while the rates across the rest of the bonds in 
the lattice are fixed to be 1 in both ways: 

•■■(-2,0) 4 (-1,0) % C (0,0) (1,0) ^ (2,0) ••• 
i ill 

It is again easy to see that the steady state density 0(f) 
now satisfies 



V 2 0(f) = -e0(O) 2<Lq 



&r , — e i 



(5) 



In the electrostatic analogue, the rhs of ([5]) corresponds 
to two oppositely oriented adjacent dipoles on the x axis, 
constituting a quadrupole charge configuration ( — h +— ) 
in the continuum limit. Using the two-dimensional 
Coulomb potential and the superposition principle, it is 
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easy to see that the density profile at large distance r 
now decays as 



(f) 



£0(0) 
2tt 



(0) 



with 0(0) = p/(l — e/2). Consequently, the particle 
current density (or equivalently the electric field of the 
quadruple) j(f) decays as r~ 3 for large r. 

In the case of an arbitrary configuration of n driving 
bonds, one uses the superposition principle to express the 
steady state profile in terms of the dipole strengths of 
the n dipoles. Using the exact expression for the Green's 
function on square lattice j§| a set of linear equations 
determining these strengths are obtained, which can be 
readily solved. 

The existence of biased bonds does not necessarily im- 
ply a breakdown of detailed balance. Localized config- 
urations of biased bonds may preserve detailed balance 
with respect to a localized potential V(r). For example, 
consider the case where all incoming links to the site (0, 0) 
are with rates (1 — e) each and the rates across the rest of 
the links in the lattice are fixed to be 1. It is easy to verify 
that the rates satisfy detailed balance with respect to the 
localized potential V(r) = — ln(l — e) Consequently, 
the steady state density has the Gibbs-Boltzmann form 
0(f) cx exp[— V(r)] leading to a flat density profile every- 
where except at the origin. 

It is interesting to note that in d = 2 dimension, the 
analogy to an electrostatic problem can be extended to 
introduce a magnetic field as well. In a general two- 
dimensional setting let denote the hopping rate from 
site i to its nearest neighbor site j. This link creates a 
pair of oppositely directed flux lines centered on the two 
plaquettes that share this link. The flux is perpendicular 
to the plaquettes, with, say, the flux to the left of the 
link points up (see FigJTk). The magnitude of the field 

generated by the (ij) bond is given by H = ln^^^. 
The total flux through each plaquette is given by the 
sum of the fluxes generated by each of its links. A neces- 
sary and sufficient condition for detailed balance to hold 
is the vanishing of the total magnetic field, as defined 
above, on all plaquettes. This is a direct consequence of 
the Kolmogorov criterion 0, [l(|. In that case, the den- 
sity profile outside the driven region is flat resulting in 
a vanishing electric field, as follows from the Boltzmann 
measure. However when the magnetic field is non-zero, 
the steady state is a nonequilibrium one, and the den- 
sity profile created by the driving bonds typically decays 
algebraically. On the other hand, bond configurations 
which result in vanishing electric field and non- vanishing 
magnetic field has a flat density profile. An example of 
such configuration is given in FigJTti. 

Another interesting density profile pattern emerges 
when one applies a global bias, say in the x-direction. 



Consider again noninteracting particles in d = 2, where 
a particle from any site (m, n) hops to a neighboring site 
with rate 1 in the north, west and south directions, while 
with rate (1 + p) to the eastern neighbor. Thus, p > 
denotes the global bias. In addition, there is the driving 
bond from (0,0) — > (1,0) where the hopping rate is set 
to be (1 + p, — e) with < e < 1. Proceeding as before, 
in the steady state, 0(f) is found to satisfy 



- p [0(F) - 0(f- ei)] + V 2 0(f) = -e0(O) 

where r — e\ = (m— 1, n) and V 2 is the discrete Laplacian 
as before. The solution can be expressed as 



0(f) = p + e0(O) g(r,0)-g(f,e 1 ) 



(8) 



where the Green's function g(f, fo) satisfies 

- M [g{r,r ) - g(f- ei,f )] + V 2 .g(f) = -8? ifo . (9) 

The solution can be obtained using Fourier transforma- 
tion. Setting (X, Y) = f — fo, one finds that for large X, 
Y 



9(X,Y) 



1 



y/ilTpX 



-liY 2 /AX 



X > 



(10) 



-pY 2 /4|*l 5 X<Q (n) 



These results have a nice interpretation as the solution of 
a diffusion equation where X plays the role of 'time' and 
Y the distance traveled from the origin. It can be directly 
seen from which, for large (X, Y), can be approxi- 
mated by its continuum version: —pdxg+dx9+9yg = 0. 
For large X > 0, neglecting the term d x , one indeed ob- 
tains an analogue of diffusion equation, /idxg = dy9 
with p playing the role of friction coefficient and X > 
being the 'time' variable and hence one obtains the stan- 
dard diffusive propagator in (fT0|) for X > 0. In con- 
trast, for negative X, one can no longer interpret di- 
rectly in terms of the diffusion equation in which the 
'time' is always a positive variable. However, upon mak- 
ing a change of variable X — > —X' and substituting 
g{X,Y) = e-» x 'h(X',Y) one gets \xd x >h = d x ,h + d Y h. 
This leads to the result (fTT]) for X < 0. Substituting 
this result for the Green's function in ([5]) one obtains 
a density profile <p(x,y) that is highly anisotropic. For 
example, for y — and as x — > +oo, the density decays 
algebraically to p as x -3 / 2 while in the direction opposite 
to the bias x — > — oo, the density decays exponentially. 
In the y direction, for fixed x, the density falls off rapidly 
as ~ exp [— /xy 2 /4|j;|] . 

The 2-d results obtained above for noninteracting par- 
ticles, with or without global bias, can be easily gen- 
eralized to arbitrary dimensions. Indeed, the solution 
in ^ holds for arbitrary dimensions d, except that the 
Green's function G(f, fo) depends on d. In the continuum 
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a long-ranged algebraic decay of the density profile 
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FIG. 2. The algebraic decay of the density profile away from 
the driving bond in the positive x direction. 



limit, the Couloumb potential G(r, rrj) behaves, for large 
\f — r |, as \r — r | _ ^~ 2 ^ for d > 2, as — In \r — f \ for 
d = 2 and as — \\r — ?q\ for d= 1. Hence, the dipole po- 
tential and consequently the density profile 4>(r) decays 
as r~( d_1 ) for large r in d > 2. In <i = I, the dipole poten- 
tial has a discontinuity at a; = 0, thus giving rise to a dis- 
continuous density profile: </>(x) = p — (e/2) 0(0) sgn(ai), 
in full accordance with Q. The results for quadrupoles 
and higher multi-poles can similarly be generalized to 
arbitrary dimensions. The analogy to a magnetic field 
discussed earlier is, however, restricted only to two di- 
mensions. 

We now show that most of the results derived above for 
noninteracting particles in presence of a localized drive 
carry through when the hard core interaction between 
the particles is switched on. We consider a symmetric 
exclusion process on a 2-d square lattice where each site 
can hold at most one particle. Our results are easily 
generalizable to arbitrary dimensions. From any occu- 
pied site (m, n) the particle attempts to hop to any of 
its neighboring sites with rate I and actually hops there 
provided the target site is empty. As in the noninter- 
acting case, we introduce the localized drive across the 
bond (0,0) — > (1,0) where the attempted hopping rate 
is (1 — e) with e < 1. It is useful to first associate an 
occupation variable T(r,t) with every site f: t(F, t) = 1 
if the site r is occupied at time t and is zero if it is empty 
at t. Clearly, the average density is (j>(r,t) — (r(r, t)). It 
is again easy to write down the time evolution equation 
of <j)(f, t) by counting the incoming and outgoing rates 
from each site. In this case the equation analogous to ([1]) 



(r) = -e(r(6)(l-r(e 1 )) 



6 r,6 



(12) 



where 0(0) in flU is replaced by C = (t(0)(1 - r(ei)}. 
While, unlike in the noninteracting case, we can not de- 
termine this prefactor self-consistently, the electrostatic 
analogy to a dipole (albeit with an unknown strength 
eC) still holds. Thus, at long distances, one still obtains 



-* eC ei ■ r 1 

*)=P~ o —+° - • 

27T r z r z 



(13) 



Consequently, the particle current density (equivalently 
the electric field due to the dipole) j(r), decays at large 
distances as in the noninteracting case up to an over- 
all multiplicative constant. In a similar way, one can also 
arrange the driving bonds so as to generate a quadrupole 
or multi-pole configurations of charges giving rise to an 
algebraic decay of the density profile with varying expo- 
nents depending on the charge configurations. 

A numerical evidence of the profiles in © and (JT3J) is 
shown in Fig. 2 where the difference of density from p is 
plotted against the distance from the driving bond in the 
positive x direction. The simulation is performed on a 
200 x 200 lattice with e = 1 and initial uniform density 
p = 0.6. The straight lines denote the theoretical results 
in ([3]) and (| 13[) with the value of 0(0) calculated using p 
and e, and C = (r(0)(l - r(ei))) = 0.3209 determined 
independently from the Monte Carlo simulation. 

In summary, we have demonstrated that in diffusive 
systems, both with and without inter particle exclusion 
interaction, localized drive can give rise to algebraically 
decaying density profiles at large distances. The problem 
of determining the density profile is mapped onto an elec- 
trostatic problem where each driving bond is represented 
by an electric dipole whose strength is determined self- 
consistently by the electric potential generated on the 
driving bond. The density profile of the driven system is 
then given by the electrostatic potential created by the 
charge distribution. An analogous quantity to the mag- 
netic field is also identified in two dimensions. 
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